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Abstract. Since close WR+O binaries are the result of a strong interaction of both stars in massive close binary 
I | systems, they can be used to constrain the highly uncertain mass and angular momentum budget during the major 

mass transfer phase. We explore the progenitor evolution of the three best suited WR+O binaries HD 90657, 
HD 186943 and HD 211853, which are characterized by a WR/O mass ratio of ~0.5 and periods of 6. .10 days. We 

■^J - , are doing so at three different levels of approximation: predicting the massive binary evolution through simple 

mass loss and angular momentum loss estimates, through full binary evolution models with parametrized mass 
transfer efficiency, and through binary evolution models including rotation of both components and a physical 
model which allows to compute mass and angular momentum loss from the binary system as function of time 

in . during the mass transfer process. All three methods give consistently the same answers. Our results show that, 

if these systems formed through stable mass transfer, their initial periods were smaller than their current ones, 
which implies that mass transfer has started during the core hydrogen burning phase of the initially more massive 
l' star. Furthermore, the mass transfer in all three cases must have been highly non-conservative, with on average 

1~ \ only ~10% of the transferred mass being retained by the mass receiving star. This result gives support to our 

system mass and angular momentum loss model, which predicts that, in the considered systems, about 90% of 
the overflowing matter is expelled by the rapid rotation of the mass receiver close to the fi-limit, which is reached 
through the accretion of the remaining 10%. 

b> '. 
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1. Introduction l|Heger fe LangeJ200(itlMevnet fe MaedeJ2 000). and evo- 



lutionary models of rotating stars are now available for 

The evolution of a star in a binary system can differ sig- many masses and metallicities . while the treatmen t f the 

nificantly from that of an isolated one with the same mass rotational processes in these models is not yet in a final 

and chemical composition. The physical processes that en- stage ( magn et ic d ynamo processes are just being included 

ter binary evolution are the gravitational and radiation | He ger et aLlBl iMaeder fc Mev^tl EoO.U they provide 

field from the companion, as well as the centrifugal force firgt ideag q{ what rQtation C£m reaUy do tQ & gtar Effects 

arising from the rotation of the system. But, most impor- of rotation; as important they are in single starS; can be 

tant, it is the evolution of the more massive component , , .i , f , , . 

' * much stronger m the components of close binary systems: 

that will influence dramatically the evolution of the sys- Estimates of the angular momen tum gain of the accreting 

tern. In certain evolutionary phases, mass transfer from gtar ^ magg transferring bina r i e s show that critical rota- 

one star to another can occur, changing the fundamental tion mav be reached quickly JPa,cketlll QRlt danger et alJ 

properties of both stars as well as their future evolution. gg| | Y()() „ k . l.an.e.lEooTTl, In order to investigate this, 

The rotational properties of binary components may we need binary evolution mode i s which include a detailed 

play a key role in this respect. The evolution of mas- treatment of rotation in the stellar interior, as in recent 

sive single stars can be strongly influenced by rotation single star models . However, in binaries, tidal processes 

~~z ; ~ T _ as well as angular momentum and accretion need to be 
bend offprint requests to: J. Petrovic, 

e-mail: petrovic@astro.uu.nl 
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considered at the same time. Some first such models are 
now available and are discussed below. 

Angular momentum accretion and the subsequent 
rapid rotation of the mass gainer may be essential for 
some of the most exciting cosmic phenomena, which may 
occur exclusively in binaries: Type la supernovae, the 
main producers of iron and cosmic yards ticks to measure 
the acce lerated expansion of the universe fYoon fe Langerl 
l2004albh . and gamma-ray bursts from collapsars, which 
the most recent stellar models with rota tion and mag- 
netic fields preclude to occur in single star s l|Petrovic et alJ 
120041 iHeger et all l2004t IWooslevI Eooi . For both, the 
Type la supernova progenitors and the gamma-ray burst 
progenitors, it is essential to understand how efficient the 
mass transfer process is and on which physical properties 
it depends. Further exciting astrophysical objects whose 
understanding is affected by our u nderstanding of mass 
trans fer comprise X-ray binaries (Chevalier & Ilovaiskv 



1992). 



1998) and Type lb and Ic supernovae IjPodsiadlowski et al 



How much matter can stars accrete from a binary com- 
panion? As mentioned above, non-magnetic accretion, i.e. 
accretion via a viscous disk or via ballistic impact, trans- 
ports angular momentum and can lead to a strong spin- 
up of the mass gaining star. For disk accretion, it appears 
plausible that the specific angular momentum of the ac- 
creted matter corresponds to Kepler-rotation at the stellar 
equator; this leads to a spin-up of the whole star to crit- 
ical r otation when its initial mass is increased by about 
20% llPackedl981^ . It appears possible that mass accretion 
continues in this situation, as viscous processes may trans- 
port angular momentum outward thro ugh the star, th e 
boundary layer, and the accretion disk (|Paczvnsyil99l . 
However, as the star is rotating very rapidly, its wi nd mass 
loss may dramatically increase l|Langerll997lll998l . which 
may render the mass transfer process inefficient. 

Observations of massive p ost-mass transfer b inary 
systems constrain this effect. lLanger et alJ (|2003^ and 
lLanger et alJ l|2004|) points out that there is evidence for 
both extremes occurring in massive close binaries, i.e. for 
quasi-conservative evolution as well as for highly non- 
conservative evolution. In the present study, we are in- 
terested in those binaries that contain a Wolf-Rayet and 
a main sequence O star. We have chosen to focus on three 
WN+O systems (HD 186943, HD 90657 and HD 211853) 
which have similar mass ratios (f«0.5) and orbital peri- 
ods (6. .10 days). As clearly the two stars in these systems 
must have undergone a strong interaction in the past, an 
understanding of their progenitor evolution may be the 
key to constrain the mass transfer efficiency in massive 
binaries: which fraction of the mass leaving the primary 
star is accumulated by the secondary star during a mass 
transfer event? 

Evolutionary calculations of massive close bina- 
ries were performed by various authors. General ideas 

about th e formation of WR+O binary systems wer e 
■ I i I I 11 I i 

given bv | Pa czvnski (1967), Kippcnhahn et alJ IjH 
Ivan den Heuvel fc Heisel l|l972l). IVanbeveren et alJ lit 




modelled the evolution of massive Case B binaries with 
different assumptions for mass and angular m omentum 
loss from the binary system. IVanbeverenT fl982) computed 
evolutionary models of massive close Case B binaries with 
primary masses between 20 M and 160 M Q . He concluded 
that most of the WR primaries are remnants of stars 
initially larger than 40 M Q and that the accretion effi- 
ciency in these systems sho uld be very below 0.3 in or- 
der to fit the observations. Ide Loore fc de Grevel l|l992h 
computed detailed models of massive Case B binary sys- 
tems for initial mass rati os of 0.6 and 0.9, assumin g an 
accretion efficiency of 0.5. IWellstein fc Langeil l|l999l) and 
IWellstein et"afl l|200ll) modelled massive binary systems 
mas s range 12..60M q assuming conservative evolution, 
and IWellsteinl l)200lh presented the first rotating binary 
evolution models for initial masses of «15M© and initial 
mass ratios q~l. 

While it was realized through these models that dif- 
ferent mass accretion may be needed to explain different 
observations, these efforts did not have the potential to ex- 
plore the physical reasons for non-conservative evolution. 

1. e., there is no reason to expect that the mass transfer 
efficiency remains constant during the mass transfer pro- 
cess in a given binary system, nor that its time-averaged 
value is constant for whole binary populations. 

It is not yet known which phy sical processes ca n ex- 
pel matter from a binary system. IVanbeverenl l|l99lj) pro- 
posed that if a binary component is more massive than 
w40-50 Mq it will go through an LBV phase of enhanced 
mass loss, which wi ll prevent the occurrence of RLOF. 
iDessart et al.1 (|2003|) investigated the possibility that ra- 
diation pressure from the secondary prevents the accre- 
tion. They found that even for moderate mass transfer 
rates (5TO _6 M0 yr _1 ) the wind and photon momenta 
which emerge from the accretion star can not alter the 
dynamics of the accretion stream. Here, we follow the 
suggestion that the effective mass accretion rate can be 
significantly d ecreased due to the spin-up of the mass 
receiving star llWellsteinl 120011 lLanger et al.ll2003l 120041: 
IPetrovic fc Langerll2004[) . 

The remainder of this paper is organized as follows. 
In Sect. we briefly discuss the observational data avail- 
able for WR+O binary systems. In Sect. |3 we derive esti- 
mates for the masses of both stars in WR+O systems for 
given initial masses and accretion efficiencies. In Sect. 0] 
we present the physics used to compute our detailed evo- 
lutionary models. Non-rotating binary evolution models 
with an adopted constant mass accretion efficiency are 
presented in Sect. |3J Our rotating models in which the 
mass accretion efficiency is obtained selfconsistently are 
discussed in Sect. El We briefly compare our models with 
observations in Sect. Conclusions are given in Sect. [HI 

2. Observational data 

There are about 20 observed Wolf-Rayet+O binary sys- 
te ms with known masse s of components in the catalogue 
of Ivan der Huchtl l|200ll) . We have chosen to model three 
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spectroscopic double-lined systems: HD 186943 (WN3), 
HD 90657 (WN5) and GP Cep (WN6/WCE)have, since 
they have similar mass ratios (g=M\yR,/Afo~0.5) and or- 
bital periods (6. .10 days). 

WN+O systems that also have short orbital periods 
are V444 Cyg, CX Cep, CQ Cep, HD 94546, HD 320102 
and HD 311884. V444 Cyg has period of 4.2 days and 
can be result of stable mass transfer evolution, but since 
mass ratio of this system is ~0.3 we did not include it in 
this paper. Orbital periods of CX Cep and CQ Cep are 
very short (~2 days) and these systems are probably the 
result of a contact evolution. HD 94546 and HD 320102 
are systems with very low masses of WR and O compo- 
nents (4M +9M Q and 2.3 M +4.1 M respectively) and 
HD 311884 is extremely massive WR+O binary system 
(51 Mq+60 Mq). Recently, an even more mas sive WR+O 
system has been observe d 83 M Q +82 M 1 Rauw et all 
l2004IPonauos et aljEool . 

The mass ratio of a binary system is determined 
from its radial velocity solution, with an error of 5-10%. 
However, to determine the exact value of the masses of 
the binary components, the value of the inclination of the 
system has to be known. Without knowledge of the incli- 
nation, only minimum m asses of t he co mponents can be 
determined, i.e., Msin 3 z. Ilvlassevl lll98l|) determined the 
minimu m mass for the WR star in HD 186943 to be 9- 
11 M Q . iNiemela fc Moffat! l|l982() determined the masses 
of the components of HD 90657 in the range 11-14 M Q for 
the WN4 component and 21-28 M Q for the O-type com- 
ponent. The masses of the WR components in HD 186943 
and HD 90657 given in T able [T] have been determined by 
lLamontagne et all l|l996f) on the bas is of improved values 
for the inclination of these systems. iDemers et all ( 2002) 
determined minimum mas ses of the componen t s of th e sys- 
tem GP Cep. Previously, lLamontagne et a 1 1 lll 996h sug- 
gested values of M WR =15M Q and M O =27M for this 
system. 

There is no obvious hydrogen contri bution in the 
WR, spectrum in any o f these system s (jMassevI Il98lt 
INiemela fc Moffadll982() . iMassevI (|l98l|) showed that hy- 
drogen absorption lines are fairly broad in the spectrum 
of HD 186943, equivalent to vsmi~ 250 km s~\ thus the 
O-type star is rotating much faster than synchronously. 

Beside the fact that the binary system GP Cep has a 
similar mass ratio and period as the other two systems, 
it has some very different properties as well. The spec- 
tral type of the WR component in GP Cep is a comb i- 
natio n of WN and W C CWN6/WCE IDemers et ai1l2002|) . 
Also, IMassevI l|l98l|) showed that, next to the main pe- 
riod of ^6.69 days of the binary system GP Cep, ra- 
dial velocities of absorption lines vary also with a period 
of 3.4698 days. He proposed that GP Cep is a quadru- 
ple sy stem, consisting of two pa irs of stars, WR+O and 
O+O. IPanov fc Seggewissl l|l990|) suggeste d that in both 
pairs o ne component is a WR star. However . IDemers et all 
(2002) showed that there is only one WR star in this 
quadruple system. 



3. The simple approach 

If the initial binary system is very close (an initial period is 
of the order of few days), RLOF occurs while the primary 
is still in the core hydrogen burning phase and Case A 
mass transfer takes place (fast and slow phase) . When the 
primary expands due to shell hydrogen burning, it fills its 
Roche lobe and Case AB mass transfer starts. During this 
mass transfer the primary star loses the major part of its 
hydrogen envelope. After Case AB mass transfer, the pri- 
mary is a helium core burning Wolf-Rayet star. During all 
this time, the secondary is still a main sequence star, but 
with an increased mass due to mass transfer. When the 
initial binary period is of the order of one to few weeks, 
the primary fills its Roche lobe for the first time during 
shell hydrogen burning and Case B mass transfer takes 
place. The primary loses most of its hydrogen envelope, 
becomes a WR star and the secondary is an O star with 
an increased mass. Case C mass transfer occurs when ini- 
tial period is of the order of years. The primary fills its 
Roche lobe during helium shell burning and mass trans- 
fer takes place on the dynamical time scale. This scenario 
is not likely for chosen systems, since some of the sec- 
ondary stars in WR+O systems have been observed to ro- 
tate faster than synchronously. This means that they have 
accreted some matter which increased their spin angular 
momentum. 

We constructed a simple method to quickly estimate 
the post-mass transfer parameters for a large number of 
binary systems for a given accretion efficiency (3. This al- 
lows us to narrow the space of possible initial parameters 
(primary mass, secondary mass and orbital period) that 
allows the evolution into a specific observed WR+O sys- 
tems. 

We considered binary systems with initial pri- 
mary masses Mij n =25..100 M© and secondaries masses 
M2 i m=25/1.7..100 M Q with an initial period of 3 days. 
We assumed that the primary is transferring matter to 
the secondary until it reaches the mass of its initial he- 
lium core (Eq. . 

Matter that is not accreted on the secondary leaves 
the system with the specific angular momentum which 
corresponds to th e secondary's orbital angular momentum 
llKino et a,lJl200lh . which is consistent with our approach 
for mass loss from the binary system (cf. Sect.QJ. Stellar 
wind mass loss is neglected. 

More massive initial primaries produce more massive 
WR stars (helium cores) in general, but if the star is in 
a binary system that goes through mass transfer during 
hydrogen core burning of the primary (Case A), this de- 
pends also on other parameters: 

-If the initial period is longer, mass transfer starts later in 
the primary evolution and the initial helium core of the 
primary is more massive. 

-If the initial mass ratio is further from unity, the mass 
transfer rate from the primary reaches higher values and 
the initial helium core mass is smaller. 
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Table 1. Basic parameters of selected WN+O SB2 binaries 



WR number 
HD number 
spectral type 
P (days) 



a sin i (Rq 
M sin 3 

M W R (Mq 

Mo (M ) 



(M ) 



WR21 
HD 90657 
WN5+04-6 
8.2546 ±0.0001 
0.04 ±0.03 
0.52 
37±3 

8.4 
50 ±4 
19 
37 



WR127 
HD 186943 
WN3±09.5V 
9.5550 
0.07 ±0.04 
0.47 
39 ±6 

9.3 
55 ±8 
17 
36 



WR153 6 
HD 211853 
WN6/WCE±03-6I 
6.6887 
0±0 
0.54 
>35.2 

73 
>6 
>21 



Notes: 

a: all parameter s from compilation of Ivan der Huchtl i200 if) , unless noted otherwise. 
bi lDemers et~ail (120021) . 



We initially want to restrict ourselves to systems that 
under go stable mass trans fer, i.e. avoid contact situa- 
tions. Iwlsteinern] {2001) found that the limiting ini- 
tial mass ratio for conservative Case A binary system 
is Mi i i n /M2 i in^l-55 and for conservative Case B sys- 
tems ~1.25. Since we allow non-conservative evolution, 
we consider initial mass ratio q<1.7 for Case A and g<1.4 
for Case B. The observed WR+O systems (HD 186943, 
HD 90657 and HD 211853) all have very short orbital pe- 
riods, between 6 and 10 days. Since, the net effect of the 
Case A±Case AB, or Case B is a widening of the orbit 
(if there is no contact), we have to assume that the initial 
periods need to be shorter than, or approximatively equal 
to the observed ones. We adopted a minimum initial or- 
bital period of 3 days to avoid that the primary fills its 
Roche lobe on the ZAMS. 

We estimated the minimum initial helium core masses, 
which are obtained by the earliest Case A systems, for sys- 
tems with a mass ratio of Afi i i n /M2,in=l-7 and an initial 
period of 3 days (Mi i j n <;41 M Q ) from the detailed evolu- 
tionary models shown later in this paper (Sect.|2J): 

M\VR,m = 0.24 * Mi )in ± 0.27. (1) 

In this linear approximation, we neglected the influence of 
the initial mass ratio on the initial WR mass. It is shown 
in Sect. 15.21 that this dependence becomes important only 
for initial mass ratios above 5~2. 

For Case B binaries, the initial WR mass does not 
depend on the initial period and the initial mass ratio 
of the system, since during core hydrogen burning, the 
primary evolves as a single star, without any interaction 
with the secondary. We estimated the relation between 
initial main sequence mass and initial WR mass as a linear 
fit from the C ase B binary systems with initial primaries 
Mi >in ^18M llWellstein fc Langerlll999(k 

M W R,in = 0.53 * Mi, in - 4.92. (2) 



The minimum initial period for a system to evolve 
through Case B mass transfer depends on the initial pri- 
mary mass and the mass ratio. We can estimate, based 
on the radii of the primaries at the end of the main se- 
quence evolution, what would be the initial orbital separa- 
tion necessary to avoid the primary filling its Roche lobe 
before shell hydrogen burning. From Kepler's law follows 
that the orbital separation is proportional to the mass ra- 
tio a^g -1 / 3 and the initial primary mass a^Mi 1 ^ 3 . Since 
q^l we can neglect this dependence and estimate the ini- 
tial period for which the radius of the primary at the 
end of MS is equal to its Roche radius. For this estimate 
we do not take into account stellar wind that will widen 
the orbit and decrease the masses. We conclude that the 
Case B limiting initial orbital period for 40 M is ^10 
days, for 45M Q ~15 days and for 75M Q ~30 days. Since 
the result of stable Case B mass transfer is widening of the 
orbit, it follows that (stable mass transfer) Case B binary 
systems can not be progenitors of the observed systems 
HD 186943, HD 90657 and HD 211853 whose orbital peri- 
ods are shorter than 10 days. However, WR star masses 
resulting from Case B evolution are practically the same 
as those from very late Case A evolution, which is still 
considered in our analysis. 

We calculate binary systems for early Case A (pi n =3 
days) and for late Case A (pin~Piimit)- 

The results are shown in Fig^for four different accre- 
tion efficiencies (/3=0. 0,0. 1,0. 5, 1.0 respectively) for early 
Case A evolution (p=3 days). Fig. [5]shows the results for 
early Case A systems (p=3 days) for all values of /3=0..1 
and for Case B/late Case A, also for all j3. 

We notice from Fig. ^ that when the assumed (3 is 
larger, the resulting WR+O systems lie further from the 
line defined by g=0.5. The reason is clear: if the accre- 
tion efficiency is higher, the secondary will become more 
massive while the initial mass of the WR star stays the 
same. Conservative evolution (Fig. ^1) produces WR+O 
systems that have small mass ratios, q~l/b.. 1/6. 



Petrovic et al.: WR+O progenitors 



5 



300 
160 

Hi 1 

£ ioo 

50 


200 
160 

£ ioo 

3 

50 



: q<0,5 .||" 


h ' t 0.1 






:c ' ' £=0.5 


: d ' ' = 1 

: * 


.0: 







10 15 EO 



5 10 15 SO 
Mwe[Mo] 



Fig. 1. Masses of both components of post-Case A mass 
transfer WR+O binary systems resulting from our sim- 
ple approach, for initial primary masses in the range 
25..100M© and an initial period of p ln —5 days, for four 
different assumed accretion efficiencies (3 (a-d). The solid 
line represents a mass ratio of q=M-wR/Mo=0.5. For an 
increasing (3, the O stars in WR+O systems become more 
massive and the WR/O-mass ratio decreases. 



The orbital period of WR+O systems depends on their 
initial orbital period, their initial mass ratio and on the 
parameter [3. If the initial period increases and there is 
no contact during the evolution, the orbital period in the 
WR+O stage will also increase. However, the orbital pe- 
riod of WR+O systems will be shorter if the initial mass 
ratio is larger. If the initial masses are very similar, the 
primary will become less massive than the secondary very 
early during the mass transfer, and afterward matter is 
transfered from the less to the more massive star, which 
results in a widening of the orbit. Conversely, the final pe- 
riod is shorter for a larger difference in initial masses in 
the binary system. 

We can draw the following conclusions: 
-The accretion efficiency during the major mass transfer 
phase in the progenitor evolution of the three observed 
WR+O binaries is small, i.e. /3=0..0.1, as for larger (3 the 
O stars during the WR+O phase are more massive and 
the WR/O-mass rations smaller than observed. However, 
we note that it is unlikely that the secondaries did not ac- 
crete at all (/3=0), since some O stars are found to rotate 
faster than synchronously. 

-The initial orbital period needs to be larger than ^3 days, 
to avoid contact at the beginning of hydrogen burning. 
-The initial orbital period should be larger than ~3 days, 
in order to obtain massive enough WR stars. 
-The initial orbital periods should be shorter than the ob- 
served orbital periods in the three WR+O systems, i.e. 
shorter than ~10 days. This excludes Case B mass trans- 
fer. 

-While the initial mass ratio Mi^n/M^in should not be 
too far from unity so contact is avoided, it should be close 
to the contact limit, since this leads to the shortest orbital 
periods and largest WR/O mass ratios in WR+O systems, 
as needed for the three observed systems. 



We conclude that if the considered three observed 
WR+O binary systems evolved through a stable mass 
transfer, a large amount of matter must have left the sys- 
tem. On the other hand, since some of the secondary stars 
in WR+O binaries have been observed to rotate faste r 
than synchronously (|Massevlll98ll: llTndernill et afl ll988l 
a certain amount of accretion may be required. 

Fig. |3 shows the resulting WR+O masses in Case A 
and Case B (latest Case A) for accretion efficiency (3=0. .1. 
If the primary star does not lose mass in a mass transfer 
during core hydrogen burning (pi n >Piimit), it will form a 
more massive WR star, as we already explained. There 
will be less mass to transfer from the primary to the sec- 
ondary, and for fixed (3 the corresponding O star will be- 
come less massive. However, since the observed periods of 
HD 186943, HD 90657 and HD 2II853 are shorter than 10 
days and Case A+Case AB widens the binary orbit, the 
initial orbital period should be shorter than observed, so 
we can conclude roughly that p m is between 3 and 10 days. 



4. Numerical code and physical assumptions 

We showed in Sect. |3that we can roughly estimate the 
parameters of the progenitor systems of observed WR+O 
binaries HD 186943, HD 90657 and HD 211853. However, 
detailed numerical models arc required in order to verify 
that the assumption of contact-free evolution can in fact 
be justified. And finally, we want to check whether the 
required mass and angular momentum loss can be repro- 
duced by our detailed selfconsistent approach. 

We are using a bi nary evolutio nary code which was 
originally developed bv lBraunl (J1998) on the basis of an im- 
plicit hydrodynamic stellar evolution code for single stars 
l|Langerlll99ll fl998j) . It calculates simultaneous evolution 
of the two stellar components of a binary system in a cir- 
cular orbit a nd the mass transfer within the Roche ap- 
proximation llKor)allll978h . Mass loss from the Roche lobe 
filling c omponent thr ough the first Lagrangian point is 
given bv lRitterl l)l98S|) as: 

M = M exp(i? - Ri)/H p (3) 
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Fig. 2. Masses of both components of post-Case A mass 
transfer WR+O binary systems resulting from our sim- 
ple approach, for initial primary masses in the range 
25..100M© and for early (pi n =3 days) Case A and late 
Case A respective Case B evolution. The assumed accre- 
tion efficiency is /3—0..1. The solid line represents a mass 
ratio of g=M WR ./M o =0.5. 

with Mo=pv s Q I \/e, where H p is the photospheric pressure 
scale height, p is the density, v s the velocity of sound and 
Q the effective cross-section o f the stream through the first 
Lagrangian point according to lMever fc Mever-Hofmei ster 
(119831) . 

Stellar wind mass loss for O stars on the main se - 
quence is calculated according to lKudritzki et al.l (1989). 
For hydrogen-poor stars (A s <0.4) we assume mass loss 
based on t he empirical ma s s loss rates for Wolf-Rayet stars 
derived bv lHamann et a l. (1995): 

log(M WR /M yr- 1 ) - -11.95+1.5 log L/L Q -2.85X,(4) 

Since lHamann fc Koesterkd (^98) suggested that these 
mass loss rates may be overestimated, we calculated evo- 
lutionary models with mass loss rate given by Eq. 01 mul- 
tiplied by factors 1/2, 1/3 and 1/6. 

The treatment of a convection an d a sem i conve ction 
which is applied here i s described in iLaneerl Jl99ll) and 
iBraun fc Langerl l)l995|) . Changes in chemical composi- 
tion are computed using a nuclear network including 
pp chains, the CNO-cycle, and the major helium, car- 
bon, neon and oxygen burning reactions. More details are 



given in IWellstein fc Lange"rl l|l999|) and IWellstein et alJ 
We use the OP AL Rosseland-mean opacities of 



Iglesias fc R.ogerslll996h . For all models, a metallicity of 



Z=0.02 is adopted. The abundance ratios of the isotopes 
for a given element are chosen to have the solar meteoritic 
abundance ratios according to Grevesse & Noels (1993). 
The change of the orbital period (orbital angular momen- 
tum loss) due to the mass tr ansfer and stellar wind mas s 
loss is computed according to lPodsiadlowski et al . (1992), 
with the specific angu lar momentum of the stella r wind 
material calculated by iBrookshaw fc Tavanil l)l993|) . 

The influence of the centrifugal force in the 
rotating models is implemented according to 
iKinnenhahn fc Thomas! ^70). The stellar spin vec- 
tors are assumed to be perpendicular to the orbital 
plane. Synchronization due to tidal spin- orbit couplin g 
is included with a time scale given by IZahnl dl977l) . 
Rotationally enhanced mass loss is included as follows: 

M/M(v Tot = 0) = 1/(1 - O) 5 , (5) 

where £=0.43, n= Vmt /v crit and v% lit =GM(l - T)/R with 
T=L/LEdd=nL/(ATrcGM) is Eddington factor, G is grav- 
itational constant, M is mass, R radius, k opacity, tw ro- 
tatin g velocity and w cr it critical rotational velocity ((Langerl 

EH. 

When the star approaches f2=l, the mass loss rate is 
increased according to the previous equation. However, 
mass loss also causes a spin-do wn of the star and equilib- 
rium mass loss rate ft eq results (|L angerll998() . If fl > Q eg , 
the corresponding angular momentum loss is so large that 
the star evolves away from the fi-limit. 

The transport of angular momentum through the stel- 
lar interiour is formulated as a diffusive process: 



dui 
~dt 



_d_ 

dm 



Uirr 2 pY iv ( 

1 am 



2w 



1 dlni 



2 dlnr 



(6) 



where v is the turbulent viscosity and i is the specific 
angular momentum of a shell at mass coordinate m. 

The specific angular momentum of the accreted mat- 
ter is determined by integrating the equation of motion of 
a test particle in the Roche potential in case the accretion 
stream impacts directly on the secondary s tar, and is as- 
sumed Keplerian otherwise I Wellsteml l|200l[) . Rotationally 
induced mixing processes and angular mome ntum trans- 
port t hrough stellar interior are described bv lHeger et alJ 
( 2000). Magnetic fields gen erated due to differential rota- 
tion in the stellar interior llSDruidl2002|) are not included 
here (however, see lPetrovic et al.ll2004|) . 

We calculated the evolution of the binary systems in 
detail until Case AB mass transfer starts. Then we es- 
timated the outcome of this mass transfer by assuming 
that it ends when WR star has ~5% of the hydrogen left 
at the surface. For this purpose we calculate the Kelvin- 
Helmholtz time scale of the primary: 



t KH = 2 • 10 7 M! 2 /(iii?ii)yr 



(7) 
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where Mi, L\ and R\\ are mass, luminosity and Roche 
radius (in Solar units) of the primary star at the onset 
of Case AB mass transfer. The mass transfer rate is then 
assumed as: 

M tr = (Mi - M WRiin )/i KH (8) 

where MwR,in is the mass of the WR star that has a hy- 
drogen surface abundance of 5%; all quantities are taken 
at the beginning of the mass transfer. We calculate the 
change of the orbital period orbit using constant value 
of / 3=0.1 for non-ro tating and /3=0.0 for rotating mod- 
els dWellsteirJl2nnl . Matter that is not retained by the 
secondary is assumed to leave the system with a specific 
angular momentum which c orresponds to the secondary's 
orbital angular momentum l|King et al.ll200lj) . 

5. Non-rotating models 

We concluded in Sect. El that massive O+O binaries 
can result in WR+O systems similar to observed the 
(HD 186943, HD 90657 and HD 211853) if accretion effi- 
ciency f3 is low. Since some O stars in WR+O binaries 
have been observed to rotate faster than synchronously, 
we concluded that /3>0.0 and assumed a constant value 
of (3=0.1 in our detailed evolutionary models. We already 
mentioned that the orbital periods of the observed sys- 
tems are between 6 and 10 days. Since the net effect of 
Case A+Case AB mass transfer is a widening of the orbit, 
the initial periods should be shorter than the observed 
ones, so we modelled binary systems with initial orbital 
periods of 3 and 6 days. 

We chose initial primary masses to be in the range 
41..75Mq. The masses of the secondaries are chosen so 
that the initial mass ratio (Mi i i n /M2 ) i n ) is <7~1. 7-2.0. An 
initial mass ratio of ml. 55 is estimated to be the limit- 
ing value for the occurrence of contact between the c om- 
ponents in Case A systems by IWellstein et alJ l|200lh for 
conservative mass transfer. Contact occurs when the ac- 
cretion time scale of the secondary (M2,acc/M2) is much 
longer than the thermal (Kelvin-Hclmholtz) time scale of 
the primary (M=Mi/iKR"), so the secondary expands and 
fills its Roche lobe. In our models, only 10% of matter 
lost by the primary is accreted on the secondary star, so 
it reaches hydrostatic equilibrium faster and expands less 
than in the case of larger (3. This is the reason why we 
adopted a weaker condition for contact formation and cal- 
culate models with mass ratios g«1.7..2.0. 

All modelled systems (except the ones that enter con- 
tact) go through Case A and Case AB mass transfer. 
Details of the evolution of all calculated binary systems 
are given in Table We discuss the details of the bi- 
nary evolution taking the system number 11 as an exam- 
ple. Fig. |3| shows the evolutionary tracks of the primary 
and the secondary in the HR diagram until the onset of 
Case AB mass transfer. This system begins its evolution 
with the initial parameters M ii j n =56 M Q , M2 i i n =33 M , 
Pi„=6 days. Both stars are core hydrogen burning stars 



(dashed line, Fig.[3J), but since the primary is more mas- 
sive, it evolves faster and fills its Roche lobe, so the system 
enters Case A mass transfer (solid line, Fig. ^5.6-10 6 
years after the beginning of core hydrogen burning. The 
first phase of Case A is fast process and takes place on the 
Kelvin-Hclmholtz (thermal) time scale (~3.1T0 4 years). 
The primary loses matter quickly and continuously with 
a high mass transfer rate (M™ ax ^3.1-10 -3 M yr _1 ). In 
order to retain hydrostatic equilibrium, the envelope ex- 
pands, which requires energy and causes a decrease in 
luminosity (Fig. EJ). At the same time the secondary is 
accreting matter and is expanding. Due to this, its lumi- 
nosity increases and the effective temperature decreases 
(Fig. |2J) . During fast phase of Case A mass transfer the 
primary loses ~19M Q and the secondary accretes 1/10 of 
that matter. After the fast process of mass transfer, the 
primary is still burning hydrogen in its core and is still 
expanding, so slow phase of Case A mass transfer takes 
place on a nuclear time scale (0.46T0 6 years) with a mass 
transfer rate of M tr ~10 -6 M Q yr _1 . After this, the pri- 
mary is the less massive star, with decreased hydrogen sur- 
face abundance. Stellar wind mass loss of the primary in- 
creases when its surface becomes hydrogen poor (X s <0A). 
At the end of core hydrogen burning the primary con- 
tracts (effective temperature increases) and thus RLOF 
stops (Fig. dotted line). When the primary starts shell 
hydrogen burning it expands (dash-dotted line, Fig. [^J, 
fills its Roche lobe and Case AB mass transfer starts. 

Fig. 0] and Fig. 03 show the evolution of the interior of 
the primary and the secondary until Case AB mass trans- 
fer. The primary loses huge amounts of matter during fast 
Case A mass transfer and its convective core becomes less 
than a half of its original mass. At the same time, the sec- 
ondary accretes matter from the primary and the heav- 
ier elements are being relocated by thermohalinc mixing. 
In Fig. El and Fig. we see the mass transfer rate and 
the surface abundances of hydrogen, carbon, nitrogen and 
oxygen. 

During Case AB mass transfer the primary star loses 
the major part of its hydrogen envelope. After Case AB 
mass transfer, the primary is a helium core burning star 
(WR) and the secondary is still a core hydrogen burning 
O star. The masses of the modelled WR stars are in the 
range from ~8..18.5M Q The orbital periods of the mod- 
elled WR+O systems vary from ~9.5 to ~20 days, and 
the mass ratios are between 0.33 and 0.53. 

5.1. Relation between initial and WR mass 

The initial mass of helium core of the primary in the bi- 
nary system depends on a few parameters: initial primary 
mass, initial period, initial mass ratio and stellar wind 
mass loss rate. If the primary loses matter due to the mass 
transfer or stellar wind during core hydrogen burning, it 
will form a helium core that is less massive than if there 
was no mass loss. If the initial period is very short, Case A 
mass transfer will take place very early in the evolution 



Table 2. Non-rotating WR+O progenitor models for (3—0.1. N is the number of the model, M\^ n and M2.i n are initial masses of the primary and the secondary, 
Pin is the initial orbital period and q m is the initial mass ratio of the binary system. t\ is time when Case A mass transfer starts, Atf is the duration of the 
fast phase of Case A mass transfer, Af™ ax is the maximum mass transfer rate, AM± t { and AM^f are mass loss of the primary and mass gain of the secondary 
(respectively) during fast Case A, At s is the duration of slow Case A mass transfer, AMi )S and AM2.S are mass loss of the primary and mass gain of the 
secondary (respectively) during the slow Case A, pab is the orbital period at the onset of Case AB, AM^ab is the mass loss of the primary during Case AB 
(mass gain of the secondary is 1/10 of this, see Sect.QJ, M W r,5 is the WR mass when the hydrogen surface abundance is X s =0.05, the WR mass at X s <0.01 
is given in brackets, Mq is the mass of the corresponding O star, q is the mass ratio M wr /Mq, and p is the orbital period of the WR+O system. The models 
are computed with a stellar wind mass loss of Hamann/6, except * Hamann/3, ** Hamann/2. 

c indicates a contact phase that occurs for low masses due to a mass ratio too far from unity, for high masses due to the secondary expansion during slow phase 
of Case A. 
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Fig. 3. HR diagram of the initial system Mi j i n =56M0, 
-^2,in=33 M , Pi n — 6 days. Both stars are core hydrogen 
burning (dashed line) until Case A mass transfer starts 
(solid line) . The primary is losing mass and its luminosity 
and effective temperature decrease. At the same time the 
secondary is accreting matter and expanding, becoming 
more luminous and cooler. After Case A mass transfer is 
finished, the primary is losing mass by stellar wind and 
contracting at the end of core hydrogen burning (dotted 
line). After this the primary starts with shell hydrogen 
burning and expands (dash-dotted line). 



of the primary, so the star will not have time to develop 
a larger core before it starts losing mass due to Roche 
lobe overflow. If the initial mass ratio {q=M \ t i n / M2, in) 
increases, the mass transfer rate from the primary star in- 
creases too and this results in less massive primaries that 
will evolve into less massive WR stars. 

In our models the primary starts losing mass by stel- 
lar wind as WR star when its hydrogen surface abundance 
goes below X s =0.4. However, the observed WR stars in 
HD 186943, HD 90657 and HD 211853 do not have obvi- 
ous hydrogen on the surface, so we assume that these WR 
stars are the result of Case AB mass transfer, with a hy- 
drogen surface abundance of X s «0.05. We also calculated 
the corresponding WR masses with A s <0.01. We plotted 
in Fig. Q3 the initial WR masses (A s =0.05 and X s <0.01) 
versus the initial primary (progenitor) masses. With 'star' 
symbols we indicated WR stars that originate from binary 



Fig. 4. The evolution of the internal structure of 
the 56 M Q primary during the core hydrogen burning. 
Convection is indicated with diagonal hatching and semi- 
convection with crossed hatching. The hatched area at the 
bottom indicates nuclear burning. The topmost solid line 
corresponds to the surface of the star. 

systems with an initial mass ratio of q~l.7 and an initial 
period p = 3 days (Tabled N 4, 8, 9, 13, 14). Large 'star' 
symbols represent WR stars with 5% of hydrogen at the 
surface and small symbols indicate WR stars that have a 
hydrogen surface abundance of less than 1%. 

We derive a relation between the initial primary mass 
and the initial WR mass (derived as a linear fit) for p=3 
days and gral.7, (A s =0.05): 

M WR = 0.24 * M 1M + 0.27. (9) 

We use this relation to estimate the initial parameters of 
the possible progenitors of the observed WR+O binary 
systems, as already explained in Sect. [21 In the same way, 
the relation between the initial primary mass and the ini- 
tial WR mass (A s <0.01) for the same systems is: 

M WR = 0.22 * M 1M + 0.56. (10) 

We also show in Fig.|^|thc initial WR masses (A s =0.05 
for binary systems N 5, 10, 11, 12, Tabled for an ini- 
tial mass ratio of ~1.7 and an initial orbital period of 6 
days (diamond symbols). We notice that the resulting WR 
masses are higher than the ones that come out from sys- 
tems with an initial orbital period of 3 days (see Sect. 15. 3(1 . 
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Fig. 5. The evolution of the internal structure of the 
33 M Q secondary during core hydrogen burning of the 
primary. Convection is indicated with diagonal hatching, 
semiconvection with crossed hatching and thermohalinc 
mixing with straight crossed hatching. The hatched area 
at the bottom indicates nuclear burning. The topmost 
solid line corresponds to the surface of the star. 

Different 'diamond' symbols for the initial primary 56 M© 
are for different mass loss rates (see Sect. 15.4(1 . Triangle 
symbols in Fig. [S] show the initial WR masses for con- 
stant initial primary mass, M 1 ; n =41 M©, but for different 
initial mass ratios (see Sect. 15. 211 

Note that the WR masses that are the result of 
early Case A progenitor evolution are significantly lower 
than ones that are the result of Case B evolution 
l|Wellstein fc Langerlll999|) . because of the mass transfer 
from the primary during the core hydrogen burning phase. 

5.2. Influence of the initial mass ratio on the WR mass 
and orbital period 

During the mass transfer phase, the mass transfer rate in- 
creases roughly until the masses of both components are 
equal. The maximum mass transfer rate during Case A 
increases with the increase of the initial mass ratio 
(Mi^n/M^in) and the resulting WR star is less massive. 
To analyse the influence of the initial mass ratio on the 
evolution of the binary system, we compared systems with 



Fig. 6. Upper plot: Mass transfer rate during Case A 
mass transfer in the binary system with M i: i n =56 M©, 
A^2,in=33M©, p m — 6 days. Lower plot: The hydrogen sur- 
face abundance of the primary (solid line) is decreasing 
during mass transfer and further due to stellar wind mass 
loss. The secondary (dashed line) recovered its original 
surface hydrogen abundance through thermohaline mix- 
ing. The primary starts losing mass with WR stellar wind 
mass loss when its hydrogen surface abundance falls be- 
neath X s =0A, represented by the dotted line. 



an initial primary mass of M li ; n =41M©, an initial or- 
bital period of p; n =3 days for five different initial mass 
ratios: 2.05, 2.00, 1.71, 1.52 and 1.37. (Table H N 1, 3, 
4, 6, 7, Fig. |5J). The system with gj n =2.05 enters con- 
tact during fast Case A mass transfer. The mass transfer 
rate in this case is very high (M«6T0 _2 M© yr _1 ), the 
secondary expands, fills its Roche lobe and the system 
enters a contact phase. The system with an initial mass 
ratio of qi n =2.00 loses ~21M© during the fast phase of 
Case A. The maximum mass transfer rate of this system 
is Mwl.8-10 -2 M© yr _1 . The helium surface abundance of 
the primary after this mass transfer is 65%, so the primary 
shrinks, loses mass through a WR stellar wind and there 
is no slow phase of Case A mass transfer (Ri< 9 R©). For 
the other three models q =1.71,1.52,1.37, the primaries 
lose less mass (~15,14,13M© respectively) during the fast 
phase of Case A mass transfer. The helium surface abun- 
dances in these systems after fast Case A mass transfer 
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Fig. 7. Surface abundance of carbon (solid line), nitro- 
gen (dotted line) and oxygen (dashed line) of the primary 
(upper plot) and the secondary (lower plot) in the system 
with Mi j i n =56M , M2,i n =33 M©, p in —& days. 



are ~30-35%. The primaries expand (i?~12-15R0) on a 
nuclear time scale and transfer mass to the secondaries 
(slow phase of Case A). 

We can conclude the following: First, if the initial mass 
ratio is larger, the mass transfer rate from the primary 
during fast phase Case A mass transfer is higher. Second, 
if the mass transfer rate is higher, the helium surface abun- 
dance of the star increases faster and if it reaches «58%, 
the primary starts losing mass with a higher (WR) mass 
loss rate and slow Case A mass transfer can be avoided. 

We also show in Fig. (lower plot) how the period 
changes during Case A mass transfer for binary systems 
N 3, 4, 6, 7. Roughly, when the mass is transfered from 
the more to the less massive star, the binary orbit shrinks, 
and when the mass is transfered from the less to the more 
massive star, the orbit widens. If the initial period is close 
to unity, the absolute difference between stellar masses 
is small, and more mass is transfered from the less to the 
more massive star during the evolution of the system. This 
results in a longer final period after Case A mass transfer. 
Systems with initial mass ratios of 2.00, 1.71, 1.52 and 
1.37 enter Case AB mass transfer with orbital periods of 
2.9, 3.9, 4.4 and 5.2 days respectively. However, the final 
period is also (more significantly) influenced by the stellar 



Fig. 8. Initial WR mass as a function of initial (pro- 
genitor) mass. Large and small symbols indicate WR 
stars with hydrogen surface abundance of A s =0.05 and 
A s <0.01, respectively. Systems with an initial orbital pe- 
riod of pin—i days and a mass ratio of q ~1.7 are indicated 
with star symbols, systems with an initial period of 6 days 
and q ~1.7 with diamond symbols, systems with an initial 
primary 41 M Q and initial period of 3 days with triangle 
symbols. The dashed line represents a linear fit for sys- 
tems with an initial period of 3 days and with A s =0.05, 
and the dotted line represents linear fit for systems with 
an initial period of 3 days and X s <0.01. 



wind mass loss rate and the amount of matter lost from the 
primary during Case AB mass transfer (see Section [OJ. 

5.3. Influence of the initial period on the WR mass 

Depending on the initial orbital period of a binary sys- 
tem, Case A mass transfer phase will start earlier or later 
in the evolution. If the period is larger, the primary will 
develop a larger core, before it starts transferring mass 
onto the secondary, and the resulting WR star will be 
more massive. To investigate the influence of the initial 
period, we compare binary systems 41 M Q -|-24M0 and 
56M Q +33M Q with p=3 days and p=6 days (Table N 
4, 5, 9, 10). 

If the initial orbital period increases for 3 days, a 41 Mq 
star will enter Case A mass transfer ~8-10 5 yr later and 
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Fig. 9. Mass transfer rate (upper plot) and orbital period 
(lower plot) during Case A mass transfer as a function 
of the change of the primary mass for systems with the 
initial primary Mi,i n =41MQ, initial orbital period pi n =3 
days and four different initial mass ratios: 2.00 (dotted 
line), 1.71 (solid line), 1.52 (dashed line) and 1.37 (dash- 
dotted line). (Table |1 N 3, 4, 6, 7) 

a 56 Mq star ~9T0 5 yr later. So, there are two things 
to point out: first, the more massive star (56 Mq) evolves 
faster, and second, a 3 days longer initial period postpones 
Case A mass transfer, for this star, by about 10 5 yr more 
than for a 41 Mq star. The net effect is a more significant 
increase of the convective core (i.e. initial helium core, i.e 
initial WR mass), for more massive star, due to the initial 
orbit widening. 

5.4. Influence of WR mass loss rate on masses and 
final period 

In our models, we assume that when the star has less than 
40% of hydrog en at the surfa c e, it s tarts losing mass ac- 
cording to the lHamann et alJ l|l995h WR mass loss rate, 
multiplied by factors: 1/6, 1/3 and 1/2 (TableEl N 10, 11, 
12). If the primary is losing more mass by stellar wind dur- 
ing core hydrogen burning, it will develop a less massive 
helium core. At the same time there will be less matter to 
be transfered during Roche lobe overflow, so the secondary 
will accrete less. 



Fig. 10. Mass transfer rate during Case A mass 
transfer for systems 41Mq+24M q (upper plot) and 
56M Q -|-33Mo (lower plot) and an initial orbital period 
of Pin=3 days(solid line) and pi n =6 days (dashed line). 
Case A mass transfer starts later in initially wider binary 
systems, the primary has more time to increase mass of 
its core and the initial WR star is more massive (Table 
N 4, 5, 9, 10). 



We show in Fig. ^] the influence of the stellar wind 
mass loss (Plot c) on the primary mass (Plot a) and 
mass transfer rate (Plot b) of systems with Mx i n =56M Q , 
-M2,in=33Mo, p in =3 and three different stellar wind mass 
loss rates (from A s <0.4): 1/6 (solid line), 1/3 (dotted 
line) and 1/2 (dashe d line) of the mass loss proposed by 
lHamann et al.l l)l995h . We notice that for higher mass loss 
rates, the slow phase of Case A stops earlier, due to the 
decrease of the stellar radius. The orbit is widening due 
to the stellar wind mass loss and the final period increases 
with the increasing mass loss rate. However, the orbit is 
more significantly widening during Case AB mass trans- 
fer. The more mass there is to transfer from the primary to 
the secondary during Case AB mass transfer, the larger 
the final orbital period. So, if the stellar wind removes 
most of the hydrogen envelope of the primary, there will 
be less mass to transfer during Case AB and the net effect 
of a higher mass loss rate is a shorter orbital period of the 
WR+O system. 
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Fig. 11. Primary mass (first plot), mass transfer rate (sec- 
ond plot) and stellar wind mass loss rate from the primary 
(third plot) until the onset of Case AB mass transfer for 
the system Mi j i n =56M0, M2,i n =33 M Q , pi n =6 days and 
three different stellar wind mass loss rates (from X s <0.4): 
1/6 (solid line), 1/3 (dotte d line) and 1/2 (dash ed line) of 
mass loss rate proposed bv lHamann et all l)l995|) (Tabled 
N 10, 11, 12). 

6. Rotating models 

When mass transfer in a binary system starts, the pri- 
mary loses matter through the first Lagrangian point (Li). 
This matter carries a certain angular momentum that will 
be transfered to the secondary. If there is an accretion 
disk, the angular momentum of the transfered matter is 
assumed to be Keplerian. If there is a direct impact accre- 
tion, like in our models, we calculate the angular momen- 
tum following a test particle moving through L\. This an- 
gular momentum spins up the top layers of the secondary 
star, and angular momentum is transfered further into the 
star due to rotationally induced mixing processes. Every 
time the secondary spins up to close to critical rotation it 
starts losing more mass due to the influence of centrifugal 
force (Eq. 5). High mass loss decreases the net accretion 
efficiency and also removes angular momentum from the 
secondary star. The secondary star is also spun down by 
tidal for ces that t e nd to synchronize it with the orbital 
motion. IWellsteinl ll200lh investigated these processes in 
binary systems with initial mass ratios close to unity and 



Fig. 12. HR diagram of the initial system Mi i i n =56 M Q , 
-M2,in=33M0, p in —6 days with rotation. Both stars are 
core hydrogen burning (dashed line) until Case A mass 
transfer starts (solid line). The primary is losing mass and 
its luminosity decreases. At the same time the secondary 
is accreting matter and expanding, becoming more lumi- 
nous. After Case A mass transfer is finished, the primary 
is losing mass by stellar wind and contracting at the end 
of core hydrogen burning (dotted line). After this the pri- 
mary starts shell hydrogen burning and expands (dash- 
dotted line). 



concluded that the accretion efficiency does not decrease 
significantly for Case A mass transfer, but in the Case B 
the parameter (3 can be significantly decreased by rota- 
tion. We present Case A rotating models with larger mass 
ratio g=Mi.i n /M2.i n =l-7..2 and find that accretion can be 
significantly decreased during Case A mass transfer. The 
reason is the following: if the initial mass ratio increases, 
so does the maximum mass transfer rate (Mmtr increases 
roughly until the masses in binary system are equal). If 
there is more mass transfered from the primary to the sec- 
ondary, the rotational velocity of the secondary is higher 
as well as its mass loss, which leads to a smaller accretion 
efficiency. 

We compare the evolution of non-rotating and ro- 
tating binary systems on the example Mi j i n =56M0 
-M2,m=33M0, an initial orbital period of p=6 days, and 
Hamann/3 WR mass loss stellar wind rate (Table N 



Table 3. Rotating WR+O progenitor models. N is the number of the model, Mi t - m and M% t i n are initial masses of the primary and the secondary, p- ln is the 
initial orbital period and qi n is the initial mass ratio of the binary system, tx is the time when Case A mass transfer starts, Atf is the duration of the the 
fast phase of Case A mass transfer, M t " lax is the maximum mass transfer rate, AM^f and AM2J are mass loss of the primary and mass gain of the secondary 
(respectively) during the fast Case A, At s is the duration of slow Case A mass transfer, AMi jS and AM^s are mass loss of the primary and mass gain of the 
secondary (respectively) during the slow Case A, pab is the orbital period at the onset of Case AB, AM^ab is the mass loss of the primary during Case AB 
(mass gain of the secondary is 1/10 of this, see Sect. 0J, Mwr,5 is the WR mass when the hydrogen surface abundance is X s — 0.05, Mo is the mass of the 
corresponding O star, q is the mass ratio Mwr/Mo, p is the orbital period of the WR+O system and Mwr,i is WR mass with A s /Ze0.01. The models are 
computed with a stellar wind mass loss of Hamann/6 except :* Hamann/3, ** Hamann/2. 
c indicates a contact phase 
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Fig. 13. The evolution of the internal structure of the 
rotating 56 M primary during core hydrogen burning. 
Convection is indicated with diagonal hatching and semi- 
convection with crossed hatching. The hatched area at the 
bottom indicates nuclear burning. Gray shaded areas rep- 
resent regions with rotationally induced mixing (intensity 
is indicated with different shades, the darker the colour, 
the stronger rotational mixing). The topmost solid line 
corresponds to the surface of the star. 



6). The rotating binary system is synchronized as it starts 
core hydrogen burning and it stays that way until mass 
transfer starts. The radius of the primary increases during 
the main sequence phase (from ~10 to ~25 R©, Fig. 119b - ) . 
but the rotation of the primary stays synchronized with 
the orbital period. This is why the rotational velocity of 
the primary also increases from ^100 to ~200 km s _1 
(Fig. H$t 0. The radius of the rotating primary increases 
faster than the radius of the non-rotating primary due to 
the influence of the centrifugal force. The result is that 
Case A mass transfer starts earlier for the rotating binary 
system (X c non «80%) then for the corresponding nonro- 
tating one(X C!rot w71%, Fig. [T3fl . 

When the fast phase of Case A starts, the secondary 
spins up (Fig. I2UL and stellar wind mass loss rapidly in- 
creases (Af sw -T(T 3 M yr" 1 , Fig.^). The accretion effi- 
ciency during this phase in the rotating system is /3=0.15 
(Table EJl- We see in Fig. QJJthat the orbital period af- 
ter Case A mass transfer of the rotating binary system is 
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Fig. 14. The evolution of the internal structure of the ro- 
tating 33 M© secondary during core hydrogen burning of 
the primary. Convection is indicated with diagonal hatch- 
ing, semiconvection with crossed hatching and thcrmoha- 
linc mixing with straight crossed hatching. The hatched 
area at the bottom indicates nuclear burning. Gray shaded 
areas represent regions with rotationally induced mixing 
(intensity is indicated with different shades, the darker the 
colour, the stronger rotational mixing). The topmost solid 
line corresponds to the surface of the star. 



shorter than for the non-rotating system (4.5 compared 
with 6.6 days). The orbital angular momentum of the bi- 
nary is changing due to mass transfer, mass loss from the 
system and spin-orbit coupling. The rotating binary sys- 
tem loses more angular momentum and the final orbital 
period is shorter than in the corresponding non-rotating 
system. Angular momentum loss in our syste ms is calcu- 
lated according to IPodsiadlowski et al.l l|l992f ) as already 
mentioned in Sect 01 and parameter a that determines the 
efficie ncy of angular momentum lo ss is calculated accord- 
ing to iBrookshaw fc Tavanil (^993). It increases with the 
mass ratio M2/M1 and the ratio between the secondary ra- 
dius and its Roche radius R2/ R\2- In rotating system the 
secondary accretes slightly more matter (/3=0.15) com- 
pared to /3=0.1 in non-rotating systems, so the mass ratio 
M2/M1 is larger in the rotating system. Second, the sec- 
ondary is spinning fast and its radius is larger than in the 
non-rotating case, and so is the ratio i?2/-Ri2- The result 
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Fig. 15. Upper plot: The mass transfer rate during Case A 
mass transfer in the binary systems with Mi^ m =56M.Q, 
-W2,in=33 M , Pin=6 days with (solid line) and without 
rotation (dotted line). Lower plot: Orbital period evolu- 
tion in rotating and non-rotating system. 



is that the angular momentum is more efficiently removed 
from the system in the rotating binary system. After the 
fast phase of Case A mass transfer, the two primaries, non- 
rotating and rotating, have almost the same mass ^34 M© 
and helium surface abundance y s ~44%. However, since 
the orbital periods are different, so are the radii of the 
primaries (~18R© for the rotating and ^23 R© for the 
non-rotating case). 

When the fast phase of Case A is finished, the non- 
rotating primary has still ~20% of hydrogen to burn 
(~7T0 5 yr), and the rotating primary has ~10% more 
than that (^1.2-10 6 yr). When the surface hydrogen abun- 
dance is less than 40%, the primaries start losing mass as 
WR stars, i.e., their stellar wind mass loss rate increases. 
Since the rotating primary has more time to spend on 
the main sequence, it also has more time to lose mass 
by WR stellar wind mass loss (7.2-10 5 yr compared with 
2.5T0 5 yr for non-rotating system). The result is that the 
non-rotating primary enters Case AB mass transfer as a 
~26 M© star with Y" s =0.75, while the rotating one is a 
~17M© star with y s =0.90. Clearly, the rotating primary 
has less hydrogen in its envelope, i.e. less mass to trans- 
fer to the secondary during Case AB mass transfer, and 



Fig. 16. The hydrogen surface abundance (solid line) in 
the primary in system with Mi,i n =56M©, M2 ; i n =33 M©, 
Pin =6 days is decreasing during mass transfer and further 
due to stellar wind mass loss. The secondary (dashed line) 
decreases its hydrogen surface abundance due to mass 
transfer. The dotted line indicates a hydrogen abundance 
of 0.4, where the primary starts losing mass with a WR 
stellar wind. 

the orbit widens less than in the non-rotating system. We 
can draw the conclusion that if rotation is included in our 
calculations, the initial WR mass is smaller and the or- 
bital period of the WR+O system is shorter than in the 
corresponding non-rotating system (Table |3J. 

We present in Fig. ^| the evolutionary tracks of 
the rotating primary and secondary in the HR dia- 
gram. Both stars are core hydrogen burning stars (dashed 
line, Fig. I12J) . but since the primary is more massive, it 
evolves faster and fills its Roche lobe, so the system en- 
ters Case A mass transfer (solid line, Fig. I12|l . The pri- 
mary loses matter quickly with a high mass transfer rate 
(M t ™ ax ft!3.2-10 _3 M© yr -1 ) and its luminosity decreases 
fFig.H2[>. At the same time the secondary accretes matter 
and its luminosity increases, but due to change in rota- 
tional velocity (Fig. I20t its radius and effective temper- 
ature are changing as well (Fig. I12H. Fig. 120b. b). During 
fast Case A mass transfer the primary lost ~19 M© and 
the secondary accreted 15% of that matter. After the fast 
mass transfer, the primary is still burning hydrogen in its 
core and is still expanding, so slow Case A mass transfer 
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Fig. 17. Surface abundance of carbon (solid line), nitro- 
gen (dotted line) and oxygen (dashed line) in the primary 
(upper plot) and the secondary (lower plot), in the sys- 
tem with Mi,in=56M Q , M 2 ,i n =33M , p ia =6 days. The 
secondary abundances are changed due to mass transfer 
of matter from the primary, thermohalinc mixing and ro- 
tational mixing. 



Fig. 18. The orbital angular momentum (upper plot) of 
the non-rotating (dotted line) and the rotating (solid line) 
binary systems with Mi j ; n =56 M Q , M2,i n =33M©, Pi n —6 
days, decreases rapidly due to mass loss from the system 
during fast Case A mass transfer and then further due to 
stellar wind mass loss. Spin period (lower plot) of the pri- 
mary (dashed line) and the secondary (dash-dotted line) 
in the above mentioned rotating binary system. 



takes place. After the primary starts losing mass with a 
WR stellar wind mass loss rate (X s <0.4) its radius will 
decrease and the slow phase of Case A stops (Fig. HSb). 
However, the primary continues expanding on the nuclear 
time scale (Fig. \19b ) and it fills its Roche lobe once again 
(Fig. El upper plot). At the end of core hydrogen burn- 
ing the primary contracts (effective temperature increases) 
and thus RLOF stops. This phase is presented in Fig. IT21 
with a dotted line. When hydrogen starts burning in a 
shell, the primary star expands (dash-dotted line, Fig.lT2l. 
fills its Roche lobe and Case AB mass transfer starts. 

The initial helium core masses are 18.6 M for the 
non-rotating and 14.8 M Q for the rotating primary. When 
Case AB mass transfer starts, the orbital periods are 7.9 d 
and 6.6 d for the non-rotating and the rotating system re- 
spectively (Fig. ll5| lower plot). The non-rotating primary 
loses ~7 Mq and the rotating one ~2 M during Case AB. 
When there is more mass to be transfered from the less to 
the more massive star in a binary system, the orbit widens 
more and the final orbital period is longer. 



Fig. 1131 and Fig. 1141 show the structure of the primary 
and the secondary before Case AB mass transfer. The pri- 
mary loses large amounts of matter during the fast phase 
of Case A mass transfer (~20 M Q ), and its convective core 
becomes less than half of its original mass. At the same 
time, the secondary accretes matter from the primary and 
the heavier elements are being relocated by thermohalinc 
mixing. Fig. 1161 and Fig. 1171 show surface abundances of 
the primary and the secondary. The secondary is accret- 
ing material from the primary and its surface abundances 
change due to this, but also due to thermohaline and ro- 
tational mixing. 

Fig. El shows the orbital angular momentum of the 
system and the spin periods of both components. The or- 
bital angular momentum of the system decreases rapidly 
due to mass loss from the system during fast Case A 
mass transfer, and then further due to stellar wind mass 
loss. The primary slows down rapidly during fast Case A 
and further due to stellar wind mass loss. The secondary 
spins up due to the accretion from the primary during fast 
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Fig. 19. Effective temperature (plot a), stellar radius (plot 
b), stellar wind mass loss rate (plot c) and rotational ve- 
locity (plot d) of the primary star in the non-rotating 
(dotted line) and rotating (solid line) binary system with 
Mi,i tt =56M , M 2:i n=33M , p in =6 days. 



Case A mass transfer and then slows down due to stellar 
wind mass loss. It spins up again during slow Case A mass 
transfer. 

The masses of modelled WR stars are in the range 
from ~11M to ~15.7M Q . Period of modelled WR+O 
systems vary from ^7.6 to ^12.7 days and mass ratios 
are between 0.35 and 0.46 (Table [|. 

6.1. Influence of rotation on the accretion efficiency 

We show in Table 03 average accretion efficiencies of rotat- 
ing binary systems during different mass transfer phases, 
and total average values with and without stellar wind 
mass loss from the primary included. During fast Case A 
mass transfer, the primary stars are losing matter with 
very high mass transfer rates (3..6.5-10 -3 M yr" 1 ). The 
angular momentum of surface layers in the secondary 
increases fast, they spin up to close to the critical ro- 
tation and start losing mass with high mass loss rate 
(~1O~ 3 M ). The average accretion efficiency during fast 
Case A in our models is 15-20%. Since this phase takes 
place on the thermal time scale, stellar wind mass loss 
from the primary is negligible during this phase. 



Fig. 20. Effective temperature (plot a), stellar radius (plot 
b), stellar wind mass loss rate (plot c) and rotational ve- 
locity (plot d) of the secondary star in the non-rotating 
(dotted line) and rotating (solid line) binary system with 
Mi, in =56M Q , M 2 ,in=33M , p ia =6 days. 

Table 4. Comparison of resulting WR masses and or- 
bital periods from non-rotating and rotating binary sys- 
tems with the same initial parameters. M\^ n , M2,i n are ini- 
tial primary and secondary mass, p; n is the initial orbital 
period, -Mwr.,5, -Mwr,i are WR masses at X s = 0.05 and 
X s < 0.01 respectively and p is the orbital period in the 
initial WR+O system where the hydrogen surface abun- 
dance of WR star is X s = 0.05. Systems are modelled with 
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Slow Case A mass transfer takes place on the nuclear 
time scale. The primary stars start losing their mass due 
to a WR stellar wind when their surfaces become hydro- 
gen deficient (X s <0.4). The WR stellar wind mass loss 
rates are of the order of : M ^-40~ 5 M Q yr -1 , and we have 
to take into account stellar wind mass loss of the primary 
during slow Case A. We calculate the mass loss of the 
primary only due to mass transfer and total mass loss in- 
cluding stellar wind mass loss, and the two corresponding 
average accretion efficiencies. If we calculate (3 S only for 
mass transfer, we notice that the slow Case A is almost 
a conservative process. The average mass transfer rates 
are ~1O _6 M yr _1 and the secondary stars are able to 
accrete almost everything without spinning up to critical 
rotation. 

Fig. shows how mass transfer rate, accretion rate 
and (3 change in the rotating model with 56M +33M , 
j»i n =6 days (WR mass loss Hamann/3) depending on the 
amount of matter lost by the primary. We also see in this 
figure the mass transfer rate from the primary in the non- 
rotating case. We can notice in the upper plot, what we 
previously discussed, that during most of the fast Case A 
mass transfer, the mass accretion rate of the secondary is 
about one order of magnitude lower than the mass loss 
rate of the primary. The primary loses ^19.3 M during 
the fast phase and the secondary gains ~2.9M Q , which 
means that on average ^15% of the mass has been ac- 
creted. However, the mass loss of the primary due to mass 
transfer during the slow phase is ~4.5M©, and the sec- 
ondary accretes ~3.9M which means that (3 S ~0.87. If 
we take into account stellar wind mass loss of the pri- 
mary stars, the average accretion efficiencies are lower. 
For example, the total mass loss of the primary during 
slow Case A mass transfer, including the stellar wind, 
in the previous example is ~1O.9M , which means that 
/3^~0.36. We neglected accret ion during Case AB mass 
transfer since IWellsteinl l)200lh showed that it is ineffi- 
cient, and since the primary stars in the modelled systems 
have relatively low mass hydrogen envelopes, and masses 
of secondary stars will not significantly change due to this 
mass transfer. (However, let us not forget that even the 
accretion of very small amounts of matter can be impor- 
tant for spinning up the secondary's surface layers and 
making it rotate faster than synchronously in a WR+O 
binary system.) Also, since this mass transfer takes place 
on the thermal time scale, stellar wind mass loss can be 
neglected. 

Finally, we can estimate the total mass loss from the 
binary systems including stellar wind, or only due to mass 
transfer, and calculate corresponding values of (3. In the 
binary systems we modelled, the primary stars lose be- 
tween 30 M and 45 M due to mass transfer and stellar 
wind, until they ignite helium in their core. The amount 
of lost mass increases with initial mass. At the same time 
the secondaries accrete 3..1OM . This means that in most 
cases 80. .90% of the mass lost by the primary leaves the 
binary system. On the other hand, the primary stars lose 
~2O..3OM only due to mass transfer, so the average ac- 



Table 5. Mass loss from binary systems. N is number of 
the model corresponding to Table /3f aBt is the average 
accretion efficiency of the secondary during the fast phase 
of Case A mass transfer. (3 s i ow is the accretion efficiency 
of the secondary during the slow phase of Case A mass 
transfer taking into account matter lost by the primary 
only due to the mass transfer, f3^ow * s ^ ne average accre- 
tion efficiency of the secondary during the slow phase of 
Case A mass transfer taking into account matter lost by 
the primary due to the mass transfer and stellar wind. (3 
is the average accretion efficiency of the secondary during 
the progenitor evolution of WR+O binary system taking 
into account matter lost by the primary only due to the 
mass transfer and (3^ d taking also into account stellar 
wind mass loss of the primary. 
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cretion of secondary stars in our models is between 15 and 
30%. 

7. Comparison with observations 

Our rotating models give generally similar results as our 
non-rotating models for (3 = 0.1. 

The rotating binary systems R6 (56 M +33 M , 
p=6 days, Hamann/3 WR mass loss) and Rl and R2 
(41M +2OM Q , p=6 days) agree quite well with the ob- 
served systems HD 186943 and HD 90657, as well as the 
non-rotating systems Nil and N12 (56 M +33 M , p=6 
days; WR mass loss rate Hamann/2 and Hamann/3). 
The system R6 evolves into a WR+O configuration with 
15M0+39M0 and p=8.5 days. I.e., its masses and pe- 
riod are close to those found in HD 186943 and HD 90657, 
even though its mass ratio of 0.38 is somewhat smaller 
than what is observed. Systems Rl and R2 evolve into a 
11M +24M WR+O system with a 9.8 day orbital pe- 
riod. I.e., period and mass ratio (0.46) agree well with the 
observed systems, but the stellar masses are somewhat 
smaller than observed (cf. Sect.|5J). Systems Nil and N12 
evolve into a WR+O system of 19 M +35 M with an or- 
bital period of 12. .14 days. In this case, both masses and 
the mass ratio (0.53) agree well with the observed ones, 
but the orbital period is slightly too large. I.e., although 
none of our models is a perfect match of HD 186943 or 
HD 90657 — which to find would require many more mod- 
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Fig. 21. Upper plot: Mass transfer (solid line) and ac- 
cretion rate (dotted line) of the rotating initial system 
56 M Q +33 M , p—6 days. Dashed line represent mass 
transfer rate in the corresponding non-rotating binary. 
Lower plot: accretion efficiency of the secondary taking 
into account matter lost by the primary only due to mass 
transfer. 

els, however, might not teach us very much — it is clear 
form these results that both systems can in fact be well 
explained through highly inefficient Case A mass transfer. 

The situation is more difficult with HD 211853 (GP 
Cep): neither the models with nor those without rota- 
tion reproduce it satisfactory. HD 211853 has the shortest 
period (6.7 d) and largest mass ratio (0.54) of the three 
chosen Galactic WR+O binaries. While we can not ex- 
clude that a Case A model of the kind presented here can 
reproduce this systems, especially the small period makes 
it appear more likely that this system has gone through 
a contact phase: contact would reduce the orbital angular 
momentum, and increase the mass loss fr om the system, 
i.e. r esult in a larger WR/O mass ratio ijWellstein et alJ 
2001). This reasoning is strengthened by the considera- 
tion that, in contrast to HD 186943 or HD 90657, the WR 
star in HD 211853 is of spectral type WN6/WNC. I.e., as 
this spectroscopic signature is not interpreted in terms of 
a binary nature of the WR component, but rather by as- 
suming that the WR star is in the transition phase from 
the W N to the WC stage (|Massev fc Orovelll989T: lLangerl 
Il99l|) . This implies that the WR star in HD 211853 must 



have already lost several solar masses of helium-rich mat- 
ter, which causes the orbit to widen. For example, sys- 
tem R6, which evolved into a 14.8 M Q -(-39.0 M Q WR+O 
system with p=8.53 days, evolves into a WC+O system 
after losing ~5 M Q more from the Wolf-Rayet star, which 
increases its orbital period by ~3 days. I.e., HD 211853 
might have entered the WR+O stage with an orbital pe- 
riod of about 4 days, which would put it together with 
the shortest period WR binaries like CX Cep or CQ Cep 
whose periods are 2.1 and 1.64 days respectively. 

During the evolution of WR+O binary system, the pri- 
mary loses mass due to WR stellar wind mass loss. WR 
stellar wind mass loss of the primary decreases mass ra- 
tio of the system and increases the orbital period, which 
means that, for example, WC+O binary system HD63099 
(M WR =9M Q , M =32M Q and p=14 days) could have 
evolved into present state through a WN+O binary sys- 
tem with q=0.5. 

8. Conclusions 

In an effort to constrain the progenitor evolution 
of the three WN+O binaries HD 186943, HD 90657, 
and HD 211853, we calculated the evolution of non- 
conservative Case A binary systems with primaries 
Mi j i=41..65 Mq and initial mass ratios between 1.7 and 
2 until the WN+O stage. We performed binary evolution 
calculations neglecting rotational processes in the two stel- 
lar components, and assuming a constant mass accretion 
efficiency of 10% for all three phases of the mass transfer, 
fast Case A, slow Case A, and Case AB. Those models 
could match two of the three systems reasonably well, 
while HD 211853, which has the shortest orbital period, 
the largest mass ratio, and a WN/WC Wolf-Rayet com- 
ponent, was found to be not well explained by contact-free 
evolutionary models: While models with shorter initial or- 
bital periods result in short periods during the WR+O 
stage, the initial WR mass is decreasing at the same time, 
which leads to smaller initial WR/O mass ratios. 

We then computed binary evolution models including 
the physics of rotation in both stellar components as well 
as the spin-up process of the mass gainer due to angular 
momentum accretion. In these models, the surface of the 
accreting star is continuously spun-up by accretion, while 
at the same time angular momentum is transported from 
the outer layers into the stellar interior by rotationally 
induced mixing processes. By employing a simple model 
for the mass loss of rapidly rotating luminous stars - the 
so called i7-limit, which was actually worked out to de- 
scribe the mas s loss processes in Luminous Blue Variables 
l|Langerlll997l) — accretion is drastically reduced once the 
star reaches critical rotation at its surface. The mass ac- 
cretion rate is then controlled by the time scale of internal 
angular momentum transport. 

Some first such model for Case A and early Case B 
have been computed by Wellstein ijLanger et alJ [2003, 
120041) for a primary mass of 15 M Q and a mass ratio close 
to one. The result was that rather high mass accretion 
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efficiencies (/3 ~0.7) could be obtained for initial peri- 
ods shorter than about 8 days. Here we find that, with 
the same physical assumptions although at higher system 
mass, the accretion efficie ncy drops to ab out 10% at an 
initial mass ratio of 1.7. As lWellsteinl l|200lj) computed one 
early Case A model for a 26 M Q +25 M system which gave 
/3=0.63, it is like the high initial mass ratio in our models 
which is responsible for the low accretion efficiency: Larger 
initial mass ratios lead to larger mass transfer rates and, 
as the time scale of internal angular momentum trans- 
port in the accreting star is rather unaffected, to smaller 
accretion efficiencies. 

Our rotating models — in which the accretion effi- 
ciency is no free parameter any more but is computed self- 
consistent and time-dependent — reproduce the observed 
WR+O binaries quite well, i.e. as good as our models 
without rotation physics, where the accretion efficiency is 
a free parameter. Our simplified considerations in Sect.|3 
have shown that this is unlikely attributable to the free- 
dom in the choice of the initial parameter of the binary 
system, i.e. initial masses and period — at least under the 
assumption that contact was avoided. In case of contact, 
various new parameters enter the model, similar to the 
case of common envelope evolution. And indeed, also our 
rotating models can not reproduce HD 211853 very well, 
mostly because it currently has a too short orbital period, 
which was likely even significantly shorter at the beginning 
of its WR+O stage. However, this of course only confirms 
the result of the simpler approaches that a contact-free 
approach does not work well for this system. 

In summary we can say that the system mass and an- 
gular momentum loss model used here — which is the 
first detailed approach to tackle the long-standing angu- 
lar momentum problem in mass transferring binaries - 
has passed the test of WR+O binaries. However, it still 
needs to be explored over which part of the space spanned 
by the initial binary parameters this model works well, 
and to what extent its results are sensitive to future im- 
provements in the stellar interior physics. The incl usion of 
magn etic fields generated by differential r otation iSjjruij 
2002) will be the next step in this direction IjPetrovic et al 
2004). 
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